Excitation hierarchy of the quantum sine-Gordon spin chain in strong magnetic field 
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The magnetic excitation spectrum of copper pyrimidine dinitrate, a material containing S = ^ 
antiferromagnetic chains with alternating g-tensor and the Dzyaloshinskii-Moriya interaction, and 
exhibiting a field-induced spin gap, is probed using submillimeter wave electron spin resonance spec- 
troscopy. Ten excitation modes are resolved in the low-temperature spectrum, and their frequency- 
field diagram is systematically studied in magnetic fields up to 25 T. The experimental data are 
sufficiently detailed to make a very accurate comparison with predictions based on the quantum 
sine-Gordon field theory. Signatures of three breather branches and a soliton, as well as those of 
several multi-particle excitation modes are identified. 

PACS numbers: 75.40.Gb, 76.30.-v, 75.10.Jm 
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Introduction. — Recently, low-dimensional spin sys- 
tems have received a considerable amount of attention. 
This is particularly due to their relevance to the quan- 
tum criticality problem, which appears to be one of the 
central concepts in modern solid state physics. Quantum 
fluctuations, significantly enhanced in spin systems with 
reduced dimensionality, give rise to a variety of strongly 
correlated states and make low-dimensional magnets an 
ideal ground for testing various theoretical concepts. To 
comprehend the role of quantum fluctuations in strongly 
correlated electron and spin systems with reduced dimen- 
sionality, it is important to explore its phenomenology in 
simple and well-controlled model systems. In this con- 
text, understanding the nature of the ground state and 
excitations in quantum spin chains is an important chal- 
lenge. An isotropic S = ^ Heisenberg antiferromagnetic 
(AF) chain with uniform nearest-neighbor exchange cou- 
pling represents one of the paradigm models of quantum 
magnetism. Its ground state is a spin singlet, and the 
dynamics are determined by a gapless two-particle con- 
tinuum of spin-i excitations, commonly referred to as 
spinons. A uniform external magnetic field causes a sub- 
stantial rearrangement of the excitation spectrum, mak- 
ing the soft modes incommensurate P, Q , although the 
spinon continuum remains gapless. Since the S = ^ 
AF chain is critical, even small perturbations can con- 
siderably change fundamental properties of the system. 
One of the most prominent examples is the S = ^ AF 
chain perturbed by an alternating g-tensor and/or the 
Dzyaloshinskii-Moriya (DM) interaction; this situation 
is realized experimentally in several spin chain systems 
0, 0, IE 0) 0- the presence of such interactions, ap- 
plication of a uniform external field H induces an effec- 
tive transverse staggered field hoc H, which leads to the 
opening of an energy gap A cx H"^^^. The gapped phase 
can be effectively described by the quantum sine-Gordon 
field theory 0, 0|. The excitation spectrum is rep- 
resented by solitons, antisolitons, and multiple soliton- 



antisoliton bound states called breathers. The availabil- 
ity of exact solutions for the sine-Gordon model allows 
very precise theoretical description of many observable 
properties and physical parameters of sine-Gordon mag- 
nets, including field dependence of^ excitation energies 
0, 0, E| and response functions [ll], This makes 

such systems a particularly interesting target for experi- 
mentally probing elementary excitations. So far, a field- 
induced gap has been observed in several S — ^ AF chain 
materials 0) S IE At the same time, experi- 
mental information on the magnetic field behavior of the 
spectrum is rather limited, and usually only the lowest 
one or two excitations can be identified. In the most de- 
tailed existing study 7] , a soliton and two breather modes 
were observed in CuCl2 • 2(dimethylsulfoxide) (CDC) us- 
ing inelastic neutron scattering. However, significant 
interchain interactions present in CDC lead to a pro- 
nounced deviation from the sine-Gordon model predic- 
tions. In ESR studies of copper benzoate 0,0, several 
modes were observed but only the lowest one was identi- 
fied. 

In this Letter we report a detailed study of the elemen- 
tary excitation spectrum in copper pyrimidine dinitrate 
(hereafter Cu-PM), which has been recently identified as 
a, S = ^ AF chain with a field- induced spin gap 0| , and is 
probably the best realization of the quantum sine-Gordon 
spin chain model known to date. Its intrachain exchange 
constant J = 36 K Q| is a factor of four larger than 
that in CDC ^E ^-nd about two times larger compared to 
copper benzoate 0. It makes Cu-PM an excellent object 
for the experimental study oi S = ^ AF chains in the low- 
temperature "nonperturbative" regime |l6l |. where the 
temperature T is small compared to the energy gap A. 
By employing high-resolution tunable-frequency submil- 
limeter wave electron spin resonance (ESR) spectroscopy, 
ten different ESR modes were resolved and their behav- 
ior in a wide range of magnetic fields up to gfisH ~ J 
was studied. By comparing the experimental data with 
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FIG. 1: The ESR transmission spectra in Cu-PM, taken in 
the Faraday configuration at a fi-equency of 429.3 GHz at 
T = 1.6 K (for explanations see the text). DPPH was used 
as a marker. 



the calculations based on the quantum sine-Gordon field 
theory, we were able to identify signatures of the three 
lowest breathers and of a soliton. 

Experimental— Cu-PM, [PM-Cu(N03)2(H20)2]„ 
(PM = pyrimidine) crystallizes in a monoclinic structure 
belonging to space group C2/c with four formula units 
per unit cell [5j- The lattice constants obtained from 
the single-crystal X-ray diffraction are a = 12.404 A, 
b = 11.511 A, c = 7.518 A, /3 = 115.0°. The Cu ions 
form chains (with a distance d = 5.71 A between neigh- 
boring Cu^+ ions at 10 K) running parallel to the short 
ac diagonal. The Cu ions are linked by the N-C-N moi- 
eties of pyrimidine, which constitute the intrachain mag- 
netic exchange pathway. The interchain Cu-Cu distance 
is 6.84 A. The Cu coordination is a distorted octahedron, 
built from an almost square N-O-N-0 equatorial plane 
and two oxygens in the axial positions. In this approxi- 
mately tetragonal local symmetry, the local principal axis 
of each octahedron is tilted from the ac plane by ±29.4°. 
Since this axis almost coincides with the principal axis 
of the g-tensor, the g-tensors for neighboring Cu ions are 
staggered. The exchange constant J = 36 ± 0.5 K was 
extracted from the single-crystal susceptibility and 
confirmed by magnetization measurements |l4l |. 

The excitation spectrum was studied using a high-field 
submillimeter wave ESR spectrometer 17]. Backward 
Wave Oscillators were employed as tunable sources of ra- 
diation, quasi-continuously covering the frequency range 
of 150 to 700 GHz. These radiation sources in combi- 
nation with the highly-homogeneous magnetic field pro- 
vided by a 25 T Bitter-type resistive magnet make the 
facility a powerful tool for studying magnetic excitation 
spectra in highly correlated spin systems. The magnetic 
field was applied along the c" direction, which is charac- 
terized by the maximal value of the staggered magnetiza- 
tion for Cu-PM [l^ . Our measurements at room tem- 



perature yield g = 2.24 ±0.02 as the effective (^-factor for 
the field applied along the c" direction, which is consis- 
tent with the data of Ref . . High-quality single-crystals 
of Cu-PM were probed using both Faraday and Voigt ge- 
ometry sampleholders. We found that the spectra ob- 
tained in the Faraday and Voigt configurations exhibit a 
pronounced similarity (possible reasons for that will be 
discussed below). 

Several resonance modes with different intensities were 
observed in experiments. A typical ESR transmittance 
spectrum obtained at a frequency of 429.3 GHz and at 
temperature T = 1.6 K is shown in Fig. ^ Within the 
experimental accuracy, the absorptions can be nicely fit 
using the Lorentzian formula for the line shape. The com- 
plete frequency-field diagram of magnetic excitations, 
collected at T = 1.6 K, is presented in Fig. |21 Absorp- 
tions denoted as Bl, S and Ul have maximal intensity, 
while the rest of excitations are about one order of mag- 
nitude weaker. It is worthwhile to mention here that the 
low-temperature behavior of their integrated intensity, 
carefully checked at several temperatures down to 1.6 K, 
strongly suggests the ground state nature of the observed 
excitations (details of the temperature evolution of the 
excitation spectrum in Cu-PM will be reported elsewhere 

Discussion. — The experimental frequency- field dia- 
gram has been analyzed in the framework of the quantum 
sine-Gordon field theory approach We used the 

following expression for the soliton gap As which is 
valid for a wide range of fields up to g^sH ^ J: 



A,, = J 
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Here c is the proportionality coefficient connecting the 
uniform applied field H and the effective staggered field 
h — cH, the parameter ^ = (2/{ttR'^) — l) , where R 
is the so-called compactification radius, and vp has the 
meaning of the Fermi velocity. Both R and vp are known 
exactly as functions oi H — qu rH I J from the solutions 
of the Bethe ansatz equations [lOj. The amplitude A^^ 
which is also a function of H , was recently computed 
numerically At a given field H, there are N = [1/^] 
breather branches i?„ with n = 1, . . .TV. The breather 
gaps A„ are given by the formula 



2AsSin(n</2). 



(2) 



At H—Q the first breather Bi is degenerate with the 
soliton-antisoliton doublet S, S. At finite H this degen- 
eracy is lifted, so that Bi becomes the lowest excitation 
and gives the strongest contribution into the magnitude 
of the gap observed in specific heat experiments . The 
sine-Gordon model predicts two more "heavy" breathers 
B2, -B3 to exist in the relevant frequency- field range. 
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FIG. 2: The frequency- field dependence of the ESR modes in Cu-PM ai T — 1.6 K. Symbols denote the experimental results, 
and lines correspond to contributions from specific excitations as predicted by the sine-Gordon theory (see the text). 



ESR probes the dynamical susceptibility x(9, at the 
momentum q — 0. However, in the present case the al- 
ternation of the g-tensor and DM interaction leads to a 
mixing of g = and q = n components. Moreover, due 
to the effective redefinition of spin operators which is 
performed to eliminate the DM interaction , and due 
to the fact that the Dzyaloshinskii vector is directed at 
an angle of about 58° with respect to the magnetic field 
0, the physical susceptibility Xphys is generally a mix- 
ture of effective longitudinal and transverse susceptibili- 
ties XsG calculated within the sine-Gordon model. For in- 
stance, the physical transverse susceptibility Xphys(Oj'^) 
is determined not only by xiG(0,<.j), but also gets a 
finite contribution from xiG'(^'^) ^'^^ x1Ig('''j sim- 
ilarly, the longitudinal component Xphysl^''-^) contains 
a small admixture of Xsci^^^) xIg^'^^^)- 
analysis of inelastic neutron experiments probing g « tt, 
this mixing could be safely neglected 0| since typically 
xiO,uj) ^ xl^i"^)) but it becomes important for the in- 
terpretation of ESR spectra, where even a small admix- 
ture of xi'^j ^) can have intensity comparable to that of 
x(0, uj). As a result, one should be able to observe contri- 
butions from xic '^gW as from xIg both in the Faraday 
and Voigt geometries. 

The energy structure of elementary excitations of the 
sine-Gordon model contributing to the low-energy spin 
dynamics is sketched in Fig. |21 Single-particle contri- 
butions to the longitudinal susceptibility xIg deter- 



mined by solitons concentrated around incommensurate 
wave vectors q = tt ± fco and by breathers concentrated 
around q = 0. Here the incommensurate shift fco = 27rm 
is determined by the total magnetization per spin m, ex- 
actly known as a function of the field; for small H one has 
Jvpko ~ gfJ-sH. In the transverse susceptibility x^g the 
dominating contribution comes from breathers at g = tt 
and solitons at g = ±/co. Thus, in an ESR experiment 
there should be several breather resonances at the ener- 
gies A„ as well as a single soliton resonance at 

Es ~ + {JvFko)^. (3) 

This exhausts the set of possible single-particle reso- 
nances (see Fig. 13); apart from them, there are various 
multiparticle continua contributing to the spectrum. 

It is important to mention that in Eqs. C3l-@ the only 
free fitting parameter is the coefficient c. Setting c — 
0.08 ± 0.002, we were able to achieve an excellent fit to 
the lowest observed mode Bl, which is described by the 
first breather gap Ai . Using this value of c, we calculated 
energies of other modes predicted by the sine-Gordon 
model, and obtained a reasonably good fit to the entire 
set of the experimental data. On the basis of the fit, we 
identify the observed resonances as follows: the modes 
Bl, B2 and B3 correspond to the first three breather res- 
onances at Ai, A2 and A3, respectively (blue lines in Fig. 
12} . The mode S is well fitted by the soliton resonance at 
Eg (the cyan line in Fig.|2Jl. This interpretation is sup- 
ported by the analysis of the temperature dependence of 
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FIG. 3: A schematic view of low-energy excitations contribut- 
ing to the transverse (left panel) and longitudinal (right panel) 
dynamic susceptibilities; S (S) denotes (anti)solitons and B 
labels breathers. 



the ESR spectra, which shows that the S mode contin- 
uously transforms into the lu — gH resonance when the 
temperature increases, in agreement with the theory |l6l| . 
The overall agreement between the sine-Gordon theory 
predictions and the experimentally obtained frequency- 
field dependence for those four single-particle resonances 
is very good. Note however that our value c = 0.08 differs 
from the result c = 0.11 obtained recently jl^ from the 
analysis of magnetization curves of Cu-PM; the origin of 
this discrepancy is not clear at present. 

The identification of the other six high-frequency ESR 
modes is more difficult since the theory does not predict 
any single-particle contributions in the relevant frequency 
region. One may nevertheless speculate that modes GI- 
GS lie very close to the edges of the soliton-breather con- 
tinua SBi, SB2 and SB3, respectively, denoted by black 
lines in Fig.|21 The edges of the BiBi and SS continua 
were not observed, presumably because they are close to 
the lines Bl and B2 and are masked by the correspond- 
ing quasiparticle contributions. For the remaining three 
modes U1-U3 we were not able to find any appropriate 
explanation on the basis of the sine-Gordon model, al- 
though we checked all possible two- and three-particle 
continua. This is especially puzzling in case of the Ul 
mode, which is one of the most intensive excitations (see 
Fig. We hope that our observations will stimulate 
further theoretical work in this direction. 

In summary, we have presented a detailed frequency- 
field diagram of spin excitations in Gu-PM, a material 
containing S = AF chains with alternating g-tensor 
and DM interaction and exhibiting a field-induced gap. 
The use of high-resolution submillimeter wave ESR spec- 
troscopy made possible obtaining a very precise infor- 
mation on the magnetic excitation spectrum. The field- 
induced gap was observed directly, and its high-field be- 
havior was studied. Ten ESR modes were resolved and 
their behavior in a broad field range up to gfisH ^ J 
was studied. By comparing the entire set of data with the 
theoretical predictions, we have provided experimental 
evidence for a number of excitations of the sine-Gordon 
theory, including solitons and the three lowest members 



of the breather hierarchy. At the same time, we ob- 
served at least one strong high-frequency resonance mode 
which does not fit into the sine-Gordon model descrip- 
tion. Our results can be relevant to understanding the 
quantum spin dynamics in copper benzoate and other 
S ~ ^ Hcisenberg antiferromagnetic chain systems. 
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